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I. INTRODUCTION. 



The important feature of the near-Earth space plasma is its irregularities [1, 2] in density, temper- 
ature of ions and electrons, currents and electric field. The existence of inhomogeneous flows and 
currents along and across Earth's magnetic field is a fundamental characteristic of space plasmas. 
One important element of the ionospheric flows appears to be a velocity shear. Compared to the 
homogeneous plasma case, sheared flows along the magnetic field introduce significant modifications 
to plasma stability. The presence of magnetic field-aligned plasma shear flows give rise to numerous 
instabilities of plasmas, which are predicted theoretically and observed experimentally in a broad 
frequency range with applications to ionospheric and fusion plasmas. The interest in scrutinizing the 
shear flow driven instabilities, that are still under intense investigation (see, e.g., recent papers [4]- [5] 
and references therein), lies, in part, in a number of observations which indicate correlation in place 
and time of broadband, low-frequency waves, and transverse ion energization with sheared flows 
on the boundaries of plasma structures [6]. Plasma shear flows were often encountered by Prognoz-8 
satellite [7] near the high latitude magnetopause which is the boundary between the Earth's magnetic 
field and solar wind. It was found[7] that such plasma flows are always associated with strong wave 
activity at ion cyclotron and lower hybrid frequencies ranges. High time resolution measurements 
of ion distribution by Fast Auroral Snapshot (FAST) satellite have revealed[2, 3] at an altitude of 
4000km narrow ion and electron beams with steep transverse spatial gradients in their energy. These 
beams were clearly identified as the sources of free energy to drive observed electrostatic and elec- 
tromagnetic ion cyclotron waves [8, 9] and particles energization in regions containing ion beams and 
field-ahgned currents. Electrostatic IC waves have been invoked in the explanations of the transverse 
anomalously strong heating of ions in ionosphere, that cannot be accounted for by frictional (Joule) 
heating [10, 11]. 

The field-aligned currents[12], ion beams[13] or relative streaming between ion species [14], 
have been proposed to provide the free energy for the development of IC instabilities. There arc some 
in situ ionospheric observations [15, 16] that support the development of the classical current driven 
electrostatic IC instability [17] in ionosphere. (Comprehensive review of the theoretical, numerical 
and laboratory experimental investigations of the current driven electrostatic IC instability was 
given in Ref. [18]). However, generally the levels of field-aligned current at which IC waves and 
transverse ion heating were observed in ionosphere were subcritical for the development of this 
instability, which has the lowest threshold for current density in the ionospheric plasma environment 
[12]. Kinetic effects, such as finite ion Larmor radius effects, electron Landau damping and ion 
cyclotron damping, are pronounced for ion cyclotron modes and kinetic theory is requisite for treating 
these modes. In Ref [19] the linear dispersion equation for plasma with the magnetic field aligned 
flows with velocities Voa (x) ||B||e2; of plasma components of a species {a — i for ions and a — e 
for electrons) with inhomogeneous number densities noa (x) was obtained using a kinetic approach 
in the local approximation, for which kxLn ^ 1 and k^Ly ^ 1, where L„ = [d\nno{x) /dx]~^, 
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Ly — [o?ln Vo {x) /dx]~^, and is the projection of the wave vector k normal to the flow velocity and 
along the velocity gradient. That equation may be presented in a form (see, also papers [4], [19] -[24]) 
as 

£ (k, u;) = 1 + ^ 5e^ (k, u;) = 1 + ^ (k, u) + (k, u)) =0. (1) 

cx=i,e a=i,e 

In Eq.(l) 5ea^ (k, u) is the conventional dielectric permittivity of the ct-species of the plasma without 
shear flow and 5ea^ (k, ou) is the velocity shear dependent part of the dielectric permittivity 5ea (k, cu). 
Here and in what follows k = {kx,ky,kz) — {k±,9,kz), where k^ — k±cos sin^, and 

kz is the projection of the wave vector on the magnetic fleld B. It was assumed in [4],[19]-[24] 
the equilibrium distribution function Fao to be a drifting Maxwellian (see below Eq.(26)), even 
though this may be somewhat unrealistic for the coUisionless ionosphere. The Maxwellian assumption 
permits calculationally convenient and comprehensive exploration of the IC instabilities. Eq.(l) was 
investigated for high frequency instabilities[19] with frequencies cu well above the IC frequency cUd, as 
well as for low frequency instabilities [20, 21] with frequencies below cjci- The investigation of the IC 
instabilities in magnetic fleld-aligncd plasma shear flows, that is a focus of this paper, was initiated 
in Rcf.[21], where IC instability driven by velocity shear of hot ion beam in plasma with cold ions 
and electrons was considered. Thereafter it was derived in Ref. [22] that plasma shear flow along the 
magnetic field may be unstable against the development of an IC instability of the hydrodynamic 
type with wavelength much less than ion thermal Larmor radius pi = VTi/oOd- That instability was 
developed even when the velocity of the electron-ion drift along the magnetic field was below the 
threshold for current driven IC instability. 

The comprehensive investigation of the IC instabilities of magnetic-field aligned plasma shear 
flow was undertaken in Ref. [4] . It was shown analytically that under conditions 

iReef^ (k,u;) | ~ (k,u;) | > jlmef (k,u;) | > \\msf^ (k,a;) | (2) 

shear flow along the magnetic fleld does not only modify the frequency, growth rate and the threshold 
of the known current driven IC instability, but it is a source of the development of the kinetic 
and hydrodynamic shear-flow-driven IC instabilities at the levels of fleld-aligned current which are 
subcritical for the development of the current driven IC instability. Shear flow along the magnetic fleld 
leads to the splitting of the separate IC mode, existing in the plasma without shear flow[17], into two 
IC modes with frequencies a;i^2 (k), which correspond to different kinds of the IC instabilities. In the 
limiting case of IC waves, which propagate almost across the magnetic field, the frequencies cl^i 2 (k), 
which arc the solutions of the linear dispersion equation 1 + e^^^ (k, cj) + ef'^ (k, cj) + Se (k, cj) = 0, 
may be obtained in the form uji^2 (k) = fc^Voi + nud + 5uji^2 (k), with \5uJi^2 (k) | <^ nuod-, |-2ji,2| = 
\5uji^2 (k) \/ \/2\kz\vTi 3> 1, \ze\ = \nujci — kz (Vgo — V^o) \/ \/^\kz\vTe % 1- In the limiting case of weak 
fiow shear, for which 

(na;„ - Kv^if A,^ [kip]) > ^'^ky^, (l - A,o {kip',) + r) , (3) 
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two kinetic IC instabilities are excited. Here and in what follows r = Tj/Tg with T^^i being the 
electron and ion temperatures, respectively, A|)^ is the Debye radius, A^n (^IPq) — In (^IPq) e~'^5.Pa^ 
and — (cTi/eBo) dlnnio{x) /dx. The frequency Re5a;i(k) = 5a;oi (k), and the growth rate, 
7oi = Im Scui (k) = + of the first instability are equal approximately to 

{nuci - kyVdi) Ain {k\pj) 



Re 5uj\ (k) = 8uq\ (k) 



{l-A,o{klpi)+T) ' 



(4) 



In 



{5u:oi (k))' 



i Lri. 



2 Fz k^Ti 



k V- 

[nwci - kyVdi) - -^-^5loqi (k) 



kzLOci 



(5) 



(5a;oi(k)) /tt (na;ci - A;^ (Vbe - Vbi) - /cyVde) / 2\ 
7ei = f; r\l ^j-j- exp [-z^) , 



where 



(6) 



nl = (m.^ - kyVd^f Al (klpj) - A^ky^^A,^ (klp^) (1 - A,o (klp^) + r) . (7) 

It follows from Eqs. (4)- (6), that 6uJi (k), defines the frequency and the growth rate of the current 
driven IC instability[17], modified by fiow shear [4]. This instability develops due to inverse electron 
Landau damping, when kz (Vog — Voi) + kyVde > ncUci under condition \ziQ\ ^ 1 of negligible IC 
damping of the mode uji. 

The solution 6uj2 (k) defines a new shear fiow driven branch of the IC waves, which is absent 
in shearless plasma fiows. The frequency Re Su!2 (k) and growth rate 702 = Im 6002 (k) — 7^2 + 7e2 of 
this second instability are equal to [4] 



Re Suj2 (k) = 6uJo2 (k) 



kzUJci {nujci - kyVdi) ' 



(8) 



7i2 



{Sujo2(k)f 2^2 T X n \ 



X 



\kz\VTi 



[nuJci - kyVdi) 



kzUJci {nuJci - kyVdi) 

2 ' 2„,2 



'^z^Ti 



k^LUd J {nuJci - kyVdi) 



(9) 



7e2 = 



{5uo2 (k))^ [Trinuci - k^ {Voe - Voi) - kyVde 



\kz\vTe 



exp {-zl) . 



(10) 



As it follows from Eq.(lO), this mode becomes unstable for any values of k±pi due to inverse electron 
Landau damping, when the velocity of the relative drift between ions and electrons is below the 



critical value Vq^^ [4], roughly estimated as Vq^' = Voi + {nujd — kyVde) /k^, i.e. under conditions at 
which the current driven IC instability modified by shear fiow does not develop. 
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In the case of a strong flow shear, for which the condition opposite to the above presented 
condition (3) is met, but with kyVl^/kzOJci < 0, the growth rates 71 and 72 are determined by Eqs.(5), 
(6) and (9), (10)), but with frequencies 5a;oi,o2 replaced with 

(1^^,34-11^)". (U, 

In the case of a sufficiently strong flow shear, when Vt^ < 0, the shear flow driven IC in- 
stability of the hydrodynamic (reactive) type is excited with frequency Su!(^h) and growth rate 7(j^) 
approximately equal to 

1 {nuci - kyVdi) Ain {k\p1) , , 

h V \ 

l-^klvl,A,^ {klpl) (1 + r) - i {ru^^ - kyVai)' (klp^) 

^<"' " i.^-AM^r) • 

This instability is in fact the extension of the hydrodynamic D'Angelo instability [25] onto the IC 
frequency range. 

An important, but still absent, element in the studies of IC instabilities of plasmas with 
parallel shear flow is an understanding of the processes of their nonlinear evolution and saturation. 
For the current driven IC instabilities of plasma without flow shear these studies were grounded on the 
nonlinear theory of the IC resonances broadening, which was developed in Ref . [26] . The saturation 
level of the current driven IC instability obtained in Rcfs. [26, 27] appears to agree with IC heating 
experiments [29, 18] . The presence of shear flows leads to more complicate picture of ions scattering in 
turbulent electric fields and resonance broadening saturation mechanism. The present work extends 
the earlier studies of the renormalized theory of the IC turbulence by including combined effect of 
plasma turbulence and shear flow, which consists in turbulent scattering of ions across the shear 
flow into the regions with a greater or smaller flow velocity and enhancing by this means transport 
of ions along shear flow. In Section II of this paper we derive the nonlinear dispersion equation, 
which accounts for the effect of IC resonances broadening resulted from the random motion of ions 
in the electric held of the IC turbulence in magnetic fleld-aligned plasma shear flow. In Section III 
we present the approximate qualitative analysis of the nonlinear evolution of the shear flow modifled 
current driven IC instability and shear flow driven kinetic IC instabilities, which resulted from the 
IC resonance broadening effect in the presence of shear flow. Finally, in Section IV we summarize 
our results. 



II. NONLINEAR DISPERSION EQUATION 

Our theory is based on the Vlasov-Poisson system of equations. We use leading center coordi- 
nates for ions X = x \ (vx/tUcj) sin 0, Y = y — (vx/tcJcj) cos0, where x,y,z are usual local particle 
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coordinates with z-axis directed along the magnetic field B, and where v± is velocity and (f) is gy- 

rophase angle of the gyromotion of ion. We find it suitable to use instead of z and new variables 
t 

Zi — z — J Vz (r) dr, 0i = + cUcit- With these variables the governing Vlasov equation describing 
the perturbation /j by the self-consistent electrostatic potential $ of the ion distribution function Fj, 
Fi — + /j, where F^q is the equilibrium function of the distribution of ions, takes the form 



dt rriiUJci V 9^ 9^ ) miV±\ dcp dv± dv± d(f) 



rrii dz\ dvz rUiUJci dY dX rrii v± d(f)i dv± rrii dz\ dvz 

This form of the Vlasov equation we consider as the most efficient for the deriving the renormalized 
solution for /j. Using the system of equations for characteristics for Eq. (14), 

dX dY dv± d(f)i dvz 



dt^ 



5$ e (9$ e oori 5$ e ujh <9$ e 9$ 



rriiUJci dYi rriiUJci dXi mi v± d(f)i rrii v± dv± rrii dzi 
e 9$ dFio e u^i dFio ^ e 9$ 9F,o ' 



(15) 



rriiUJci dY dX rrii 50i dv± rrii dzi dv^ 
the following nonlinear solution for the function fi with known Fjo is obtained : 



m J 



1 dFio CJ^ dFio ^ dFio 



dt'. (16) 



Ud dY dX v_L 901 dv_L dzi dv^ 

Supposing that the particle orbit disturbance 6X due to the electrostatic plasma turbulence 
is sufficiently small, we find the solution for system (15). From the first equation of that system we 
obtain 

t 



X = X + dX, dX 



rriiUJci 



f 

-JW'*- (IT) 



where X and Y arc the guiding center coordinates averaged over the turbulent pulsations. All other 
equations of system (14) have the following approximate solutions: 

t 

Y^Y + 5Y, 5Y^ / -^dh, (18) 

rriiUJci . oX 



t 

v±^v± + 5v±, 5v± = / dti, (19) 

m v± J 86 



t 

+ 50, H = -——I^dt,, (20) 

rrii v± J ov± 
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e r (9$ 



Vz = Vz + Sv:„ Sv^ = I -^dti. (21) 

m J oz 







In Eqs.(17)-(21) the perturbed electrostatic potential $ is defined in variables X, Y, v^, 0, Vz and 

5X, 5Y , 5vj^, 5(f), 5vz as 

$(r, i)= ^ / dkJn i ^^^^ j exp ik^X + ikyY + ik^z 

n=—oo \ ct / 

+ikz {vz + Vio {^)) t-in{^- Udt - 9) exp {ik.5v {t)) j duj^ (k, a;) exp {-iuut) , (22) 
where the perturbations of the ions orbits due to wave-ion interactions, 

k5r (t) = kJX {t) + ky6Y {t) + kzSz (t) - ^^^'"^ sin {(f) - 9) - cos (0 - 9) 50 {t) , (23) 



with 



t 

Sz 



(t) = J 5vz (h) dh + K (^) y" (h) dh, (24) 

are included. In Eq.(24) the first term corresponds to the scattering of ions along the magnetic field 
by the Ez — —ikz^ projection of the turbulent electric field. The second term corresponds to the 
combined effect of turbulence and shear fiow. This effect consists in turbulent scattering of ions 
across the magnetic field along the velocity shear from the fiow layer with velocity Vq (X) into the 
layer with velocity Vq [X) + Vq (X) 5X and enhancing by this means the transport of ions along the 
shear fiow. The order of value of the ratio of the second term in Eq.(24) to the first one is kyV^ j kzOJd- 
In laboratory and ionospheric plasmas[6] this parameter may exceed unity considerably. In Eq.(24) 
the terms of the second order in 5X, 5Y , 5v^, 5vz and S(f) are omitted. The Fourier transformed 
Poisson's equation, 

$ (k, a;) = ^ ejria (k, a;) = ^ e„ ^ f^dv, (25) 

a=i,e oc=i,e 

and Eq.(16) in which the potential $ (r, t) is related with $ (k, co) by Eq.(22) with 6r determined by 
Eqs.(17)-(21), compose the system of nonhnear integral equations for $ (k, lo) and fa- We obtain from 
this system the nonlinear dispersion equation which accounted for the scattering of ions by turbulence 
in shear fiow. We assume the equilibrium distribution function F^q to be a drifting Maxwellian 



4i 4i 



^ (O 2^3/2 -2 I ■ (26) 
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Then the ion density perturbation will be equal to 



Srii 



oo oo 



(k, ^) = y dy^fi k, v±, (f), v^) 

^ oo oo 

2m^ ^ j dvz j dv±v±J rfa;i$(k,a;i) J dte^^™)*^ dtiFio{X,v) 

' "=-°°-oo -oo 

X exp ^iuji {t - ti) - ik^ [v^ + Vio (X)) {t - ti) - inuj^i {t - h) 



X exp 



-ikAr {t,ti) 



2 f k±v± \ 



1 7 — 7 — 

kyVdi H kyVz TlLOci f^z'^z 



OJr. 



(27) 



where Ar [t, ti) — 5r (t) — 5r (ti). As in conventional renormahzed theory [26, 27] we account for the 
average effect of the perturbations of ions orbits in Eq.(15). For this we can use the cummulant 
expansion [2 8], 



(28) 



exp (-ikAr {t, ti)) ) = exp f ^ — ^ - ikAr {t, ti) 



where {■■■)c is the cummulant. In the following we use the simplified approximation that the particles 
scattering by plasma turbulence is a Gaussian process, for which expansion (28) under Markovian 
approximation is reduced to a single term [2 6], 



exp (—ikAr {t,ti)) ^ ~ exp ^~^( ~ ^i) 

With these approximations the ion density perturbation will be equal to 

drii (k, ^) = y dvfi {uj, k, v±, 0, Vz) 



(29) 



{X) J2 {klpl) {ky'^di - nuci) $ (k, a;) Ri (k. 



a; 



' (^) E (^i^') kyT^ f 1 - ) ^ ^) R2 (k, U) 



kzOJci 



(30) 



where 



oo 

Ri {k,(jj) = J rfrexp i(u - nud - kzVio (X) - ^^klv^^T^ - ^((k(5r (r) )' 



(31) 



oo 

i?2(k,a;) = J drrexp i(u - nujc - kzVio (X)^t - ^klv^iT^ - ^(^{k5r (r))' 



(32) 



with r = t — ti instead of ti, are resonance functions which replace the familiar resonant denominators 
of the linear theory. Using Eq.(23) and omitting the terms which oscillate with frequency nujci 
{n — ±1, ±2, ...) we can write the mean square displacement as 

I ((k ■ 6rf) = U{k- 6rf)^ + k,kz {6X6z) + kykz {6Y6z) + {{6zf) , (33) 



where 



<(k • 6rf)^ = kl {{6Xf) + kl {i6Yf) + 2k,ky {6X6Y) + 



1 t2,,2 1 1.2 



2 a;2 



The term ((k • (^r)^)g exists in plasma even without shear flow[26]. The first three terms in Eq.(33) 
describe the diffusion of the ion guiding center coordinates X and Y and the last two terms describe 
the random changes of the IC radius v±/ujci and phase angle <p of the Larmor orbit. In the asymptotic 
limit of large time compared with the correlation time Tcorr of the turbulent electric field along the 
particle orbit, and after averaging over ions velocity assuming a Maxwellian distribution we obtain 
that the term ((k • ^r)^)^^ is equal to 

- ((k • 5vf)^ = C^t = :r^Re J2 / (kl) l^e-'^i^? [2 (fc.fci, - k^^kyf [kip]) 



^2 ^2 



oo 

/ 



X / cirexp 



-i - nuc - k^ViQ {X))t- -^{{^i^r (r) )' 



(35) 



Last three terms in Eq.(33), which are due to scattering of ions along the magnetic field, are equal 
to 



with 



k^k, {6X6z) + kyk, {6Y6z) + hi {{Szf) = Cst^ + C^t^ 



Co 



Re 



oo „ 



ni=— oo 
oo 



dkl|$ (kl) fe-'^-'^'k^, ( 1 + ) {k^k^y - k,,ky) 



xln^ {k\pl) J drexp 



{uj - nujc - k^Vio {X))t- - (^ (ki5r (r) ) ' 



(36) 



(37) 



3m- 



Re ^ / dki|$ (kl) \^e-'''^fhlkl (l + 

ni=— oo ^ 



xL 



ni 



oo 



drexp 



1 



■i {u; - nwc - k^Vio (X) ) r - - ^ (ki^r (r) ) ' 



(38) 



where the term with C2 originates from the first two terms on the left of Eq.(36), whereas the term 

t 

with C3 resulted from the last left term of Eq.(36). Due to the term Vq (X) / 5X (ti) dti in Eq.(24), 
term C3 includes the effect of the anomalous viscosity, resulted from the turbulent redistribution of 
the shear flow moment. 

The frequencies and growth rates of ion cyclotron instabilities, which were considered in 
Ref.[4], are determined under condition of weak IC damping, i.e. for |a; — /c^V^o (X) — ncUcil/kzVTi 3> 1. 
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Under this condition the asymptotics of integrals (31) and (32) are determined by the contributions 
over the time interval < r <^ |a; — /c^V^o {X) — nujci\~^ and contributions from the stationary phase 
point. Our calculations of the integrals Ri and R2 give that the contribution from the stationary 
phase point is negligible for both integrals and main contributions proceed from the vicinity of r = 
point and are equal to 

00 

R,^ [dr exp USujT - C^r' - ^st^) ^ f + ^ - , (39) 



00 

R,^ [ drre.p (i5u;r - C.r^ - C,r') ^ --^ - 7^ - (^0) 

J ^ ' [buj) [buj] [buj] 

where now ^uj (k) is the renormalized version of the linear ^lo (k), presented in Introduction, 

?>u (k) = u; - k^ViQ {X) - nujd + iCi, (41) 

and C2 = C2 — \k1v'^^. In Eq.(39) the first term i/5uj includes known [26] ion cyclotron resonance 
broadening effect; the following two ones, which are proportional to C2 and C3, reflect the combined 
effect of turbulent scattering and shear flow. Using above obtained results for Ri and R2 in Eq.(30) we 
obtain the renormalized version of the equation for perturbed ion density, bni (k, cj), which accounted 
for the the combined effect of shear flow and ions scattering by electrostatic turbulence. 



5ni (k, uj) = (k, uj 



1- E 



UJ 



- k^VQ [X) - kyVdi + iCi 



n=—oo 



5uj 

(1 + An - ^f^* (k, 1^ E 

$ (k, u) 2^ —-3^ An {k^Pi ) ( C2 + 3z— 1 (42) 

n=— 00 V / \ ' 



n=— 00 



Inserting Eq.(42) to Poisson's equation (25), we obtain the nonlinear dispersion equation 

1 + £oi (k, U) + £e (k, u) + Esh (k, uS) = 0, (43) 

which together with Eqs.(35)-(38) comprises the system of nonlinear integral equations for spectral 
intensity |$(ki)p and frequency uj. This system is derived under condition that the dominant 
nonlinear effect is a randomization of ion orbits. This effect enters the dispersion equation only 
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through resonance functions (39) and(40). In Eq.(43) we introduce the notations 

£oi (k, u) — Re e\^^ (k, ou) + Reef^ (k, u) 

1 k V ^ 2 

" ^ E Trf^^iniklpl), (44) 



/c2A2j. /C^^ci ^fr'oo (^^) 

1 



where z^o — {uj — /c^Veo) /v^^zfTe and 



+ A£e(k,a;), (45) 



s..(k,a;) = ^ E --^ A4^.P.)(^^2 + 3^^j 

The system of equations (35), (37), (38) and (43) is the main result of this theory. To obtain the 
solution oj (k) of Eq.(43) we expand it about the solution ujq (k) of the linear dispersion equation, 
l + £oi(k,u;o (k)) + 1/A;2A2,, = 0, 

1 + £oi (k, (k)) + (a; + iCi - (k)) + As^ (k, a;) + Ssh (k, a;o (k)) = 0. (47) 

The solution is 

(k) = C.0 (k) - iC, - A^e(k,^)+ e..(k,a;o(k)) _ 

agp^ (k, a;o) 

Thus, the nonlinear growth/damping rate, which accounts for the scattering of ions by plasma 
turbulence in along-field shear flow is equal to 

7 (k) = Im w (k) = 70 (k) - Ci + ^sh (k) , (49) 

where 

ImAee(k,cJo(k)) Ime^fe (k, cjq (k)) 

^»« = — fec»(k,a,o) ■ ^-w^-^^ik;^- (5") 

In the case of the shearlcss plasma flow IC resonance broadening effect, determined by Ci, leads to 
the saturation of the current driven IC instability [26]. Below we examine the role of the shear flow 
driven term '-jgh (k) in the nonlinear evolution of the shear flow modifled and shear flow driven kinetic 
IC instabilities. It follows from Eq.(50) that when (k) < the combined effect of the shear flow 
and IC resonance broadening leads to the steady state, 7 (k) = 0, and the saturation level of the IC 
turbulence is estimated from the equation 

7o(k) = Ci-7,^(k). (51) 

However, when 7s/j (k) > 0, a nonlinear instability may occur. 
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III. NONLINEAR STAGE AND SATURATION OF THE 



ION CYCLOTRON INSTABILITIES OF MAGNETIC FIELD- 
ALIGNED PLASMA SHEAR FLOW 

In this section we apply the renormahzed dispersion equation (43) to the quahtative analysis of 
the nonlinear evolution of the kinetic IC instabilities presented in Introduction, which are devel- 
oped due to inverse electron Landau damping. It follows from Eqs.(6) and (10), that growth 
rates 70(1,2) — 7e(i,2) of the investigated instabilities, attain theirs maximal values at \ze\ — 
\nujci — kz {Veo — Vio) \/V^\kz\vTe 1- For nwd ~ \kz (T4o — Vio) \ that gives the estimate 

kz VTi 1 , . 
, (52j 

k±_ VTe k±Pi 

which links possible values of kz/k± and k±pi for which the maximum values of the growth rates 
may be developed. Applying this estimate to condition (3) of weak flow shear, and to the condition 
of strong flow shear, opposite to (3), we come to the following estimates: 



for weak flow shear, and 



< 



OJr. 



> 



VTi k±Pi (1 - AiQ {k\p1) + r) 

A-in (^IPi ) 



Vt, 



(53) 



(54) 



VTi k±pi (1 - Aio {k'ipl) + r) 

for strong flow shear, which determine the possible values for kj^pi in the cases of weak or strong flow 
shear for which maximal values of the growth rates may be attained. Accounting for these estimates, 
as well as the results of the linear theory presented in Introduction, we consider the cases of plasma 
flows with weak and strong flow shear separately. 



A. Weak flow shear 

In this subsection we consider the case of weak flow shear, limited by condition (3), for which two 
kinetic IC instabilities may develop [4]. The flrst one is the shear flow modifled current driven IC 
instability, which is excited when kz (Voe — VQi)-\-kyVde > nuJd with frequency ujqi — nuJci — kzV(yi-\-5ujQi 
and with growth rate 701 = 7ii + 7ie, where 5ujq\, and 7ie are determined by Eqs.(4), (5), (6), 
respectively. The second one is the shear flow driven IC instability [4]. It is excited with frequency 
ujq2 — nuici — kzVoi + SuJo2 and with growth rate 702 = 72i + 72e; where Su!o2, 721 and 72e are determined 
by Eqs.(8), (9), (10), respectively, when the velocity of the relative drift between ions and electrons 
is below the critical value Vq^\ roughly estimated as Vq^'' = Voi+ {nujd — kyVde) /kz, i.e. under 
conditions at which modified by shear fiow current driven ion cyclotron instability does not develop. 
In the case of shearless current the dominant nonlinear effect responsible for the saturation of the 
current driven IC instability is the turbulent broadening of the ion cyclotron resonances [2 6, 27], 
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determined by the term Ci in Eq.(49). The shear flow introduces new combined effect of shear flow 
and turbulent scattering, which is incorporated in ^ygh term in Eq.(49). Here we evaluate the joint 
effect of shear flow and IC turbulence, determined by the terms Ci and in (49) on the nonlinear 
evolution of the above mentioned shear flow modifled and shear flow driven kinetic IC instabilities. 



1. Shear flow modified ion cyclotron current driven instability 

Under the condition of the weak flow shear (3), the ratio \joiC2/Cs\ of the C2 and C3 contained 
terms in jsh is 



701C2 



y 



It follows from Eq.(55) that for k^pi > 1 I701C2/C3I > A;^p| > 1 and for k^pi < 1 I701C2/C3I > 
{kyPi^"^^ » 1. For this reason, for /c_LPi » 1 as well as for /c^pj <g 1, term with C3 in Eq.(45) may 
be omitted and is determined by C2 term in the flrst line of Eq.(46) as 

Note, that because of h^VTi <^ ^c^oi, the second line in Eq.(46) may be neglected. Now we estimate 
the ratio Cxj^ahi^) under condition of weak flow shear. We have that 



for k±pi 1 as well as for k±pi <^ 1. Therefore the effect of shear flow on the saturation of the 
shear flow modifled IC current driven instability is sub dominant. The level of the IC turbulence in 
the steady state, e$/Te, where ^ — (/ |$ (ki) pdki) is the root- mean-square (rms) magnitude of 
the perturbed electrostatic potential, is estimated from the balance equation 70 (k) = Ci. Using the 
mean value theorem for the integral over ki in Eq.(35) we have 

>f^^V" (68) 

in the short wavelength, k±pi 3> 1, part of the spectrum, and 

~ n-l 

— - 2 2 {k^pi) > 2 2 ^59^ 

-tj yVTe^ci J 

in the long wavelength part, k^pi <^ 1, where we have used Eqs.(52) and (53) in limits k^pi ^ 1 
and k^pi <^ 1, respectively. It follows from estimates (58) and (59) that $ has a very low level (58) 
for IC turbulence in the short wavelength part of spectrum and attains its maximal value for n = 1 
and k^pi < 1 [26] , which is equal to 

e$ 1 
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2. Shear flow driven ion cyclotron instability 

Now we estimate the ratio 702C2/C3 which determines relative importance of the C2 and C3 contained 
terms in Imegh (k, cc;o2 (k)). It follows from Eqs.(7), (9), (37) and (38) that 



702C2 



~ r 



kz LOci 



Arn{klpi)kl (l-Ao(^ip?)+^ + ^^Al,J' ^ > 



where constraint (3) was used. Therefore for long wavelength, k_i_pi < 1, IC mode ijJ2 (k) the term 



with C3 is dominant in IvuSsh- In this case 7^^'' is equal to 



2 rwjci + 24- 



where the approximation of the separate IC mode was used. In Eq.(62) C3 is determined by the 
relation 

which follows from (38) and (39). In Eq.(63) the first term on the left side, C3, is dominant over the 
second one when 

S ^ ky^, kyVl^ 1 

% - ul kz.,.Al\k\pl) ^ ^ 

However at level (64) 7^^-* ~ 24 {ky^/^QijkzLOcif' ^{i2 ^ 702, i-e. the saturation occurs before the second 
term in left hand side of Eq.(63) becomes comparable with the first one. In the case k^pi < 1 



(2) T \ h V . 

Thus the nonlinear damping rate 7^^^ (62), which for k±pi < 1 may be approximated as 

(2) ^ 24{r + klpj)C, 
^^'^ ~ 1 + r A;44,' ^ > 

dominates over Ci in this case. The level of the IC turbulence in this case will be determined by 
the balance equation 7^^'' ~ 702 with C3 term determined by Eq.(63) without the second term. This 
balance equation gives the level 

where Eq.(3) was used. Level (67) gives the estimate for the amplitudes of the perturbed potential 
for different values of the wave number k for long wavelength IC waves. This level is determined by 
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balancing linear growth rate with nonlinear damping rate, which we consider in our estimations as 
formed by perturbations of the electrostatic potential with spectrum width Ak ~ k about k. The 
ultimate saturation level at which a steady state occurs over all wave number space is determined 
as a level which is sufficient for the stabilization of the wavenumber region of the modes with the 
maximal growth rate. It follows from Eqs.(lO) that the growth rate 702 — 7e2 attains its maximal 
value at Z2e 1- Thus the ultimate saturation level is determined by Eq.(67) with kz/k± determined 
by Eq.(52) and and it is equal to 



VTiUJci (fcjopf ) 

In Eq.(52) k±oPi is determined by Eq.(53) and it is equal to 



1) 



(68) 



V VTe<^ci 



l/(2n-l) 



(69) 



Level (68) is in maximum for the fundamental IC mode with n = 1, for which 

1 > ko±p^ > (70) 

The condition \zin\ > 1 of the weak ion cyclotron damping, under which Eq.(42) was obtained, 
imposes other restriction on the admissible values for ko±pi, which for mode 5u!2 is 



1 > ko±pi > 



(71) 



The maximal level of IC turbulence, for which conditions (70) and (71) are met both, is of the order 
of 



1\ 



(72) 



for 



and 



< 



1) 



VTi 
VTe 



1/2 



VTi 
VTe 



(73) 



(74) 



when the condition opposite to (73) met. Note, that levels (72) and (74) appear to be much lower 
then level (60) for modified by shear flow current driven IC instability. 

Now consider the short wavelength, k±pi > 1, part of the spectrum of IC waves. Accounting 
for Eqs.(52) and (53) we obtain that 



702C2 



kz t^r.i 



An (klpj) kl - VTeVk 



> 1. 



(75) 
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which furnishes the dominance of the C2 terms in ^gh for that part of the spectrum. Nonhnear 
damping rate ^sh is equal now to 

6702C2 , . 

Ish ~ jz 72- (76) 

The estimation of the ratio Ci/^sh with ^sh determined by (76) gives, that 

— ~ -— ~ — — > 1. (77) 

-fsh k±pi kz VTi yU 



Thus, the term ^{^h in Eq.(49) may be omitted. The ampUtude of the electrostatic potential $ in the 
steady state may be determined in this case from the balance equation 702 (k) = Ci, which yields 
the estimate for the rms amplitudes of the perturbed potential for different values of k^^pi for short 
wavelength IC waves. 



(78) 



Ti VTeOJci{k^Pif'^' 

where the Eq.(52) was used. Eq.(53) and condition of weak IC damping of IC waves impose on k^pi 
the restriction 

1/2 

^^y, J >^.P.>^ (79) 

For that range of k^pi level (78) is in the range 



VTe \i^cij ^ Ti ^ \VTe'^ciJ 

and the electric field strength E±_ — —k±^ is 

E,<^(^}^'\ (81) 

As a demonstration of level (81) we take the data detected in magnetopause by satellite Prognoz-8 
[7, 22], 5 = 2 ■ lO^-'^G, VoJiOci = 0.5, ^ Ti = lOOeV, and obtain the estimate E± < 0.4 ■ lO'^V/m 
at fundamental cyclotron mode, which is in good agreement with measured value[22] E± ~ 0.5 • 
lO-^V/m. 



B. Strong flow shear. Shear flow driven ion cyclotron instabilities 

Under conditions of sufficiently strong flow shear, for which condition opposite to (3) is met, two shear 
flow driven IC kinetic instabilities may be developed for kyVQ/kzUJci < with frequency determined 
by Eq.(ll). These instabilities are excited with a growth rate determined by (5) and (6) when 
kziy^e — Vqi) > nUci, and with growth rate (9) and (10), when /c^ (Voe — Voi) < nuJci- For short 
wavelength, k±pi > 1, ion cyclotron waves the estimate 

1/2 

> 1 (82) 



701,02^*2 



rik^Pif^ 
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is valid, which defines the dominance of the C2 terms in ^gh for that part of the spectrum. Nonhnear 
damping rate ^sh in this case is equal to 



(84) 





ky Vq^ 







The estimation of the ratio Cxj^sh with 7^/1 determined by (83) gives that 



:i+t) 

24t 



(A;±Pi) 



1/2 



ky UJci 

kz Voi 



> 1. 



Therefore the term jgh in Eq.(49) may be omitted. The amplitude of the electrostatic potential $ in 
the steady state may be determined in this case from the balance equation 70 (k) = Ci, which gives 
the following estimate: 



1\ 



1/2 



(85) 



k±Pi \k±u, 

By using Eqs.(52) and (54) we obtain from Eq.(85) the ultimate estimate for the value of the per 
turbed potential in the short wavelength part of the spectrum, 



(86) 

VTei^ci/ 

Now we consider the long wavelength, k±pi < 1, spectrum subrange of IC waves. Comparing 
the contributions of C2 and C3 to 7^/1 under condition of strong flow shear we obtain that 



701,02C2 kzVTi k] 



OJri kl 



kz ^ci 
ky 



1/2 



(An (klp^))'^' (1 - Ao (klp^) + r + eXl,) 



1/2 



< rklp^ (1 - (fclp?) + r + k^Xl,) 



1/2 



(87) 



Thus for long wavelength, k±pi < 1, ion cyclotron waves the nonlinear growth rate 75^ in the case of 
strong flow shear is dominated by the C3 term and is approximated as 



24 



7,, (k) ^ Y^CsA^ (klp^) 



ky^Oi 



kzLOci 



(88) 



where the coefficient C3 is determined from the equation 

.2 



C3 + 



2e' 



m: 
2 



[dk,\^{k,)\'A,^{kl^p^)klkl( 



klyViO 



e 



ni=—oo 
00 



klzUJci J {SUJ+) 



, ^ J rfk,i$ (ko \^A^ {ki^pi) kiki (^i^y 



701 



ni=— 00 



(89) 



Remarkably, 7s/i (k) determined by Eq.(88) is positive and it may be responsible for the development 
of the nonlinear instability. From the estimate of the ratio \Ci/'^sh\ under condition of strong flow 
shear (54), 



Ci 



Ish 



r^(l + r) 
24 



k.UJr. 



kyVoi 



< 



r^(l + r) kiv. 



2„,2 

Ti 



24 



«1, 



(90) 
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it follows, that the dominant nonlinear term in Eq.(49) in the long wavelength part of the spectrum of 
IC waves in the case of strong flow shear is 7^^, determined by Eq.(88) with coefficient C3, determined 
by Eq.(89). Thus the combined effect of strong shear flow and turbulent scattering of ions, determined 
by the ^sh, leads to the nonlinear instability with growth rate (88). Note, that the first term of the 
left-hand side of Eq.(89) dominates over the second when 

f<rA!^ (91) 
and it becomes less than the second term at the level of turbulence above (91). At the level 

^ > r^'^Jf (92) 

the nonlinear growth rate becomes greater than the linear growth rate 70 (k). 



IV.DISCUSSION AND CONCLUSION. 

In this study, we have presented the renormalized theory of the IC turbulence in the magnetic 
field-aligned plasma shear flow. The developed theory extends the earlier studies [26, 27] of the 
renormalized theory of the IC turbulence by including a new combined effect of plasma turbulence 
and shear flow, determined in nonlinear dispersion equation (43) by the term Ssh (k, a;). This effect 
consists in turbulent scattering of ions by the IC turbulence across the shear flow into the regions 
with a greater or smaller flow velocity and is the manifestation of the anomalous viscosity due to ion 
cyclotron turbulence. Analytically it manifests in nonlinear broadening of IC resonances. We have 
derived the approximate solution (49) of that equation and use it for qualitative analysis of the role 
of the discovered effect in nonlinear evolution of the shear flow modified current driven IC instability 
and shear flow driven kinetic IC instabilities. We have considered the limits of weak and strong flow 
shear for short wavelength and long wavelength parts of the spectrum of IC waves separately and 
arrived at the following conclusions. 

1. The shear flow modified current driven IC instability, which develops under condition of 
weak flow shear (3), saturates as the ordinary IC instabihty driven by shearless current on the high 
level (60) in the long wavelength, k±pi <^ 1, part of the spectrum and on a very low level (58) in the 
short wavelength, kj_pi ^ 1, part. The effect of the shear flow, determined by the term jsh (k), is 
negligible on the saturation of this instability. 

2. The long wavelength part of the spectrum of the shear flow driven IC instability under 
condition of weak flow shear saturates due to combined effect of shear flow and turbulent scattering 
of ions, determined by the term ^fsh (k) in Eq.(48), on levels (72) or (74) (depending on the condi- 
tion (73)), which appear to be much lower then the corresponding level (60) for current driven IC 
instability modified by shear fiow. 
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3. The saturation of the short wavelength spectrum subrange of the shear flow driven instabihty 
under conditions of weak (49) , as weU as strong (49) flow shear arises from the turbulent scattering 
of ions by IC turbulence as in shearless plasma and determined by term Ci in Eq.(49). It occurs at 
very low levels (78) and (86), respectively, which are comparable with corresponding level (58) for 
shear flow modifled IC current driven instability. 

4. Nonhnear evolution of the long wavelength part of the IC turbulence spectrum developed 
by the shear flow driven instability under conditions of strong flow shear is determined by the ^sh (k) 
term in Eq.(49). In this case combined effect of shear flow and turbulent scattering of ions introduces 
a principally new effect into the nonlinear development of the IC turbulence, which is absent in 
shearless plasma flows. It is the shear flow driven nonlinear instability with growth rate (88), which 
above level (91) becomes greater than the corresponding linear growth rate 702 (k), determined by 
Eqs.(9) and (10). This effect resembles the effect of the negative viscosity. It contrasts dramatically 
with effect of ion scattering by IC turbulence in shearless plasma flows[26, 27], where it leads only 
to saturation of IC instability. 

It has to be noted here that we consider the model of the collisionless plasma. However even 
weak ion-neutral collisions, roughly 2-10^^ of the IC frequency, are stabilizing. Ion collision frequency 
may be absorbed into the linear growth rate, leading to reducing the maximal growth rate. When ion- 
ion or ion-neutral collision frequencies arc even marginally greater than an ion cyclotron instability 
growth rate, such instability does not develop [31]. The effect of ion collisions on the evolution 
of the considered ion cyclotron instabilities in near-Earth space plasmas determines the threshold 
altitude for the instabilities development. Our collision-free analysis is completely applicable to shear 
flows formed by solar wind around the Earth magnetosphere, to auroral acceleration region where at 
an altitude of 4000km the NASA's FAST (Fast Auroral SnapshoT) satellite observed quite narrow 
ion beams[2] with Vq/uJci ratio of the order of [6] 0.3 for hydrogen ions and of the order of 10 for 
oxygen ions. The collision-free analysis of our paper is pertinent to the F-region of the ionosphere for 
altitudes above 180km. In this region ion-neutral coUision frequency Vin is roughly 4 • 10~^ of the ion 
cyclotron frequency and this value diminishes with altitude. That value may become less than the 
linear growth rate 7e2 (Eq.(lO)), 7e2 ~ (^^T^/^^Te)^ (K/^ci)^ (^IPi Mm (^IPi )) ^ci, and is much less 
then 7± ~ {vTi/vref^'^ {Vl/udf^ {k±piAin {k^pD^^Ud (see Eqs.(lO) and (11)) for r ~ 1, k^pi ~ 1 
for shear flows with Vl^ujd- 

This paper is the flrst attempt of the analysis of the nonlinear evolution of the IC turbulence 
of magnetic fleld-ahgned plasma shear flow. Here we have considered only the effect of scattering of 
ions by the IC turbulence in shear flow on the nonlinear evolution of kinetic IC instabilities. It is well 
known that other nonlinear effects, such as quasilinear flattening of the electron distribution function 
in wave-particle resonance region[17], and weak turbulence effects[32] may change the ultimate picture 
of the nonlinear stage of the IC turbulence. Real process of the nonhnear evolution of the instabilities 
considered will include the interplay of all these processes. That defines the need in the development 
of additional nonlinear theories for the plasma shear flows, which, however, are absent today. These 
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are IC weak turbulence theory for plasma shear flow, which has to include the processes of the induced 
scattering of IC waves on ions and electrons and processes of the decay of IC waves, ordinary and 
renormalized quasilinear theories for ion and electron distribution functions. Adequate coverage of 
such a broad area of nonlinear plasma physics is difficult in a limited-page article and much more 
have to be done for the development of a comprehensive nonlinear theory of the turbulent state of 
plasma shear flows. 
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